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ABSTRACT 


Properties of a stationary sequence of random variables ee 
which have exponential marginal distributions and random linear 
combinations of order one of an i.i.d. exponential sequence te.) 
were discussed by Lawrance and Lewis (1976); they called this model 
the EMA] (exponential moving average of order one) point process. 
This paper will investigate the estimators of the parameter 8 of 
the EMAl process, and some basic properties of the EMA2 process, 


and then extend these results to the EMAKkK process. 
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tee LNTRODUCTION 


Properties of the stationary sequence of positive random variables 
i} which are formed from an independent and identically distributed 


exponential sequence fe} according to the linear model 


BE Veep sooability 6B, 


Pn eats Ee 
2Be Peal with probability 1-8. 


were discussed by Lawrance and Lewis [Ref. 1]. They gave a fairly com- 
plete picture of this model, and called it the EMAl1 (Exponential Moving 
Average of order 1) process. It is clear that the adjacent elements of 
this sequence are correlated, but that the dependence is no greater than 


order one, 1.e. x. is independent of Xx, en, orate «ana .SO LOren: fOr 


+2 1+3 


pairs and triples. 

In this paper, methods of estimating 8 and the properties of the 
estimates of § will be discussed, and then the properties of an 
analogous second order process are investigated. The new process, called 
the EMA2 model, is a sequence of positive random variables tx, } defined 
by 

bao), es 


2 


Po = + eo Ul Hie Ale —1 0) penal ee a 
4 Son eee W.p ( BR, ae is ’ r ) 


+ + -p. (l1- 1- : ale} 
an eee e wep. ( oe. 8) ( 


oe 


The purpose in the creation of this model is to provide models for 


data with longer dependencies than that obtained with the first-order 





model and to examine any tendencies of the upper bound on the serial 
correlations to increase. For the EMA] model O<p,<1/4 and P,.=0 £Or 
fae see eee Or the EMAZ2 model it is shown that 0<p,,0,S1/4 and 
Pp, =0 Robe K—Ipis. Li Fact, the {x. } form a sequence of exponential 
random variables, and it will be seen from (1.1) that the successive 


2: ey aaa a will be correlated. This model is also an 
i¢+l th ee 


elements Xin xX 
alternative model to a renewal process. 

The EMA2 model is shown to be a stationary point process. Distri- 
bution of the sums of x. are discussed, and the joint distributions 
of two adjacent intervals x. are derived and appear to be new bivari- 
ate exponential distributions. Extensions of the model and estimation 
problems are briefly discussed. 


In developing the properties of the process, the similarities to 


a backward second order moving average which is defined as 


mre: wW.p. 8 ; 
Ze 2 
mee-(f8 €,+8 €, Veo ee (0-8 6 Slei=0,21,....); 
ai es a es ecm | gt Dat al 
? eee a 5. (1k) (1-8 ) . Cie 
ad eae ae. Eo 2 1 


will also be pointed out. Properties of the processes are very similar, 


but those of the forward model (1.1) have simpler derivations. 





Pies enter REVinW OF THE EMAL PROCESS 


The EMAl model iS a stationary point process with exponential mar- 
Samal Gistribution of the intervals 1X. 4 - Further Xs is dependent 
on X. 1 and Xa , but independent of all others, so the correla- 


a O = GO (x ae 
fon Pp, =c¢ ( i 


441? = o) (ES Is) 7 P,.=0 1 ee eae ca route eter 


The Laplace transform of the p.d.f. of THX tk tee FX, is 





r ee 2Bs) iL 
#2] el, 2a) 
A+S8 ((A4+BS) [A+ (14+B8) Ss] 
Let Noo be the number of events occurring in the interval (0,t] be- 


Sapming at an arbitrary event; and let F(t) denote the distribution 


or T +; then 
r 


(on a > 
ProbiN, =cJ=F (t) F 41 (t)- a6 


with P(t) =1 pomeeeOy., 6ne p.d.f. Of Non gives the generating 


Gimeclon as 


n (£) ay 


E[z © )-= = 5 (23 wom s TTF (t)-F_ 4 (t)}=1+(z-1) Ez F(t). (2.2) 
r=0 r+ r=1 a 


Miserting (2.1) in the Laplace transform of (2.2) gives 


x 8 (1+8)s7+[-B(1-B) z+2B+1]AS+A? 


Y (258) =—<—<— $< (2.3) 


£ (sth) [8 (148) S24 (1428-282) As+(1-z) A7] 
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Peererenelating (2.3) with respect to z , then setting z=1 , gives 


the Laplace transform of the intensity function m,(t) 


, as 
x feCr + ps) (A+ (i+8)s] 
<2.) CoS 
3 B (1+8)s (A+s) [stA/(87+8) ] 
and inverting (2.4) gives 
2 iat 
188) 6 mats, (6 ONES At, (52 aeas 


m (t)={ 
f r(1483Ate >t) (B24+8=1) 


ine, JOint distribution of Xx. and Xa is a bivariate exponen- 


tial. Using a double Laplace transform we get 


£ 
X. X47 'Sy0 S4)=¥ (Bs) [By (Bs) +(1-B) ¥(s, +85, )) [8+ (1-8) p(s.) ] 


A*(A+Bs, +85.) 
is (A+B8s,) (Ats,) (Ats, +Bs,) " (2.5) 


and using triple Laplace transform gives 


KKK 


f. 


Be 9 fae ( 


1 (84752153) 
Bee) 164(0,)+(2-999 0, +80,)1 199 (80,) +(1-819(0,480,)9 [8eL-£0¥(0,)1 
Mitferentiating (2.5) with respect to S. 7 SeLetang S,=0,, and 


inverting with respect to s and then dividing by the marginal (exponen- 


tial) density of Xe , gives 


a eee ache ee Gla b) t/P 
E(X,|X,_j=t)=) ~ [BAt+ i ie 


11 





Similarly, 


~(1-B) AE/By 


-1 
ogee) =A) (1th-e 


The two conditional variances are given by 


fee Zeno me zest isAe) _-(1-B)xt/B_ __8? o72 (1B) At /B} 


Pome x, =t)=X “| ea ce eae 
ayy 72 pi+8+62-83 5, _B | At, _-h(1-8)t/B__-2(1-B)At/B 
Paix, ,=t)=A ~ { = 21 Se ee ae e Me 
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III. ESTIMATING 8 IN THE EMA] MODEL 


The EMA1 model iS not time-reversible, and this comes out clearly 
in higher order joint moments. The results lead to a method for esti- 


mating 8 in the EMA] model. 


Berine 2 2) 
Cc rot (x) . FEN Gl ate : 
1,2 REC, Xe EO ECKL), 


2 Z 
a Oe.) X 
1 SRIFECK XL) E(X )E(X. 0), 


which when k=l gives 
2 


Cg Q)HE (KX )-E(X. Pa) ; 


2. : 
Coy CLFE(X X, 1) E(X{)E(X, 3): 


By the construction of EMAI1, we have: 


5 Bree W.p. Bae 
ae = (C262 t=O, 212 ,.-..) 
Be epgece e,  1ee e? fal WD. (1-8) 4. 


Hence, uSing straightforward combination, we get the joint expectation 


of = and Xa as 
B(X¢X, ,, )=E( p3c2e, |) B45 Boer ey te €40545)8 (8) 
oe 
+E (BP eFe; 1428763 CF ar tB OSes ot 2BE Es Pa 40 +Be3 tet yy 549) 8)”. 


Simplification of this result leads to 


CS )= 3 (2448- -28%-283) which implies that S , {= cules sh) (Aas) 


1S 





Similarly, we get 


28 


Z = = 
, (1) eo . BY (148). 


L a 7] ) , 


1, 


Therefore, if we let 


C , 2) 


> : 


we have a function of 8 which decreases monotonically from 2 when 
fee to 3/2 when 61 . Thus there is a unique solution for 8 for 
emyediven r ; note that when $8 is 0 or 1, the ratio is not defined. 


Solving (3.1) we get 


c i: 
each he 





C by, = 
er; IDA 2c, 4) CF a? 
FSS aa a 1)-c (1) 
ean | Coy! Ve 
1 
C1 2! ) 
For estimating § , define 
nel 
2 1 eo -. -2 
<5, 2s ain! ae a) } 
1=1] 
n-1 
a j 2 -2. ,7 
Sa ee 2 Xia 
1=1 ' 
ae a 2c, ,(1)-c, , ) 
Cy, )-c) 4) 


Now we check all the estimators, to see if they are asymptotically 


Maolased or not. 


=] , 1 n n 
LBM XG, Else Dx, Dx), 
=) 1=1 1i=1 


a 1 n 
er 


uk 
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Examining the estimate of the product of the means we have 


n n 
2 ia 2 
* = Tt a a eet SEP os 
} X ) X= (Xj 4X54. x) CX tX, ) 


=X, 42 (n-1) (X.%s FX Xi 47 )+(n-1) (n ito. 
Tis 
n n “re 
ial 2 if i 2 3 2 
cep = 6nt+2 (n-1) (4468 -28°-48")4+2n -6n+4]) 
Boe EST ;? Pears or ( ) 
1-1 j= A 
2 4 D 3 
m2 Spf w= -(2438~-B2-28 2) + —z 
Siac ( B-B e3 


--2 
ana when no , E(XX y32/r3 
. . 2 eZ - 
Thus, since the estimators of E(X.X. 14) and of E(X,X. 11) are unbiased, 
we get that asymptotically, 


N 


2 : 
Ec, ,()I= “ys 8-8) +8); 


. “Le 
Pe aa Sed le yx 8-8) (248). 


1.e. Both of these are unbiased estimators when n is large. 
geeewe now look at the ratio estimator of 8 , namely 8B to see if it 


is asymptotically unbiased. Note that the denominator in the expression 


“nN 


for 6 is identically zero if B=0 or f#=1 , so that in what follows 
we assume that O<f<l . let 


ray A 


2c, 5 (1) -C., 1 )=¥ ange let 


, , 


Se) oy 9) -2- Then 


ea. an 1 Dee Cen: 
pt =enic, 51) E[c, , (J a (28 26> ake and 


wee KIC eee ea 
i Bic, | Qy-Bic, 0) I= —5(28-287)-n, 


HO y/u OW mb. 
Now we can write (Y/Z)= oS ry (ie cia I where Y'=Y/yu_, so that 
yu 2/1 Ul y 
u 4 zZ 2 u Zz 
ae 7— 2 
B(¥/Z)= 7 ice) 148(Y"(— Meee ob | ya 


Z Zz Zz 


we assume that conditions for this expansion to hold as m© are met. 


iS 





Since ae and ety’ )=l1 , i1£ B=Y/Z is to be unbiased, we 
must show the rest of the terms in (3.1) are all zeros. 


Thus look at 
1 (Zu) /u=E (YZ) Juul. 


we have 


an an “A ™. 

n-l n-l 5 n-l n-1 
ria ct> 2 X. Pan Reeecpopene YexexX. (> 2X X51! 
i=l i=l Rell Seal 

n-l 2 n-l 5 

Me ee 2 
=[ (2U-w) (W-U)]/(n-1)*-[ 2X, 2 Xf (W-U)I/[n*(n-1)], 
7=]1 ~“i=l] 


where = = =X.X* 
U de Bee 2 aly Wee U. x eq! 


=) eu .= Lee TmCueyr =X =X. 2 
gud oe ek ay ile Aig Ee al 


In addition (2U-W) (W-U) =3UW-2U7-WwW?. 


Further we get 
n-1 n-l n-l n-2 


NS mee ed: U We ee le ae Ue 
1 gag Pi-2 ly a i-V yly itd 


7 UE W. a. iMin3- 


Oe ce UU (m= 3) (n-4) U5 Us, 3, 


+2 2 W.W. +2 % W-W. +(n-3) (n- 4)W. We 
; + 
i=l 70 er | Bae 1 it2 oes 
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It can be shown that all these joint expectations have finite expected 
value so, when m© , those terms only with coefficients n will go to 


: ae 
zero when multiplying by aap Mee thus asymptotically, 


E[ (2U-W) (W-U) / (n-1) 2] =2 (3UW-2U7-W~) / (n-1) 2 ae x 
3E(U;)E(W;43)-2E (Ui) E (Ui, 3) -E(W;) E (Wy, 4)=(48-4874+48°-88%4+485) /)§, 


ae) 
i+1 


E (W;) =E(X£X, 4) =(2+48-287-289) /25, 


since E(U;)=E(X,X =(22R=2G-)/X°, 


and 
i oa ee, 3 
Similarly, E[ © X. 2 X¢(W-U)/(n~-n) ]-------- > 
isl Sia 


2 
E(X,)E(X.) (E(W,)-E(U,)]=(48-487) /a°, 


3 
Mmenece, E(YZ)= == (48 -88°+48°) , when n is large. 
A 
Asymptotically, 
E[Y' (2-1) /u,1=(48°-88"+4B°) /(u ju A°)-1=0. 
Za Ne ae, 
In (3.1), the rest of the terms in the braces will also approach zero 
when nis large, so E(f8)=B—, i.e. 8 is an unbiased estimator. 


An alternative way to estimate 8 is to use 
a A 
B=[A°C | (1)1/29 2 


where 6 is an estimator of 1 the first order serial correlation 


au 
of EMA1, and p,=8 (1-8) Reem, 1. ps5] 


Define a n-1 _ 2 
oes XSi ag / (ne) - G0) 
Then me : n-l n 5 
E(p )=A ee PEN Ey) }/n7}. 


an 





Again using the same argument as above, we have 





= m2 2 
E(X.X,,,)=(1+B-B°) fA". 
Also n > Nn 9 ia) NS sole al 
exe) =. 2 Xs+2 ¥ XX. t2 2 2; re Saree 
i=l + i=] ied ei? eS 
so that E[ (x) “J=en72[2n/d242 (n-1) (146-82) /22 
no = 
+(n-1) (n-2) /X*]-------- fay 
x 5 noo 
ga -auentily E(p )=A*{ (n-1) (1+8-8*)/[A* (n-1) ]-E[ (x) [ae eass eA oil Oleas)) ae 
Thus 2, is an unbiased estimator for O, when n is large. 
Now assume n° and let Y=)"C, 5 (1) and Z=2, Tas 
xe ZB 4S BXI+B); and E(Z)=28(1-8), 
aes 4 (Be -8°-8°4+8°), 
and by the expansion used above 
u fai 4-\ 
ez) =—"{e (y') -B[v' (——*) ]+E[¥'(——4) 2]-....}3 (3.2) 
Me Ue ae 


We want to show that the terms in (3.2) beyond the first are zero, so 


we look at 
E(¥' (2-v,)/u,J=E (YZ) /pyyz-l. 


We have n-l] n 
2G) ae ee _ 445 Z 
1=1 =] 
n-1 n n 
Z 
[ 2 El, (0-1) = a, ae» x~/n2). 
= 7 fei) She & 


Ife. 





Therefore 


5 2 2 
EiG.2) ----- aX [E(XX4 J ECKGX, HE (XX, E(X, E(x) 


“E(X,)E (Xj) B(X;X2 43) +E (X,) (Xj) B (Xz) B84) 
=4(87-6°-B*+8°), 


so that, 


; -]=1-1=0. 
Ety! (2-u_) /u J SE (2) Ny by ies lieaal 


Similarly, we can show the rest of the terms in the braces of (3.2) all 


approach zero when n is large. Hence 
A 
E(B) =E(¥/Z)-l=1+8-1=8, 


ra 


aA 
meeeweo is also an unbiased estimator of 8 whenn is large. 


iS 





IV. COMPARISON OF THE ESTIMATORS 


“~~ 
van 


It has been shown in the last section that the two estimators B, B 
are unbiased asymptotically, provided that 840, or B41, but when the 
sample size n is not large enough, the bias term should be considered. 
For simplification, any finite term divided by the second or higher power 
of n will be neglected. 


(1). All the esti- 


For 8 plete Y=2C (1) -C,, bi Z=C (1)-c 
td 


lie area 


Mators here are defined as before. Thus 


ee 28 148) /.°-4 (243B~B*~28°) /nd°= au; 


B(Z)=28(1-8)/A*=p,; 
ae 
E(Y/Z)=n uy {1-{E (YZ) /u su -1) +... 03 


Using those results listed in APPENDIX B and neglecting the higher power 


terms, we have 


E(YZ)=E(Y¥)E(Z) + (-12-56084+1472B7-11488 ?+5328"%-11208° 


+7408%+24887-1648%) Ben . 


Hence E(B) =(2u u,-E (YZ) ]/us 
1 1 
= B +} « ——-~ (124+5288-1488874+12128°-516BR" 
487(1-8)* n 


My 


+1088B°-7408°-2488'+4+164B°%). (4.1) 


a 3A A 
F = bd — ie 
Seep), tet Y=A oT ae Z=20, 
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Again all the estimators here are defined as before. Thus 


B(Y)=28 (1-8) (1+B)-4 (2+38-B7-2B8°)/n=ny and 


E(Z) =28 (1-6) + (2+48-48*) /n=u,- (4.2) 


2 
In (4.2), the maximum value of (2+4B-48') occurs at B=1/2 , and 
equals to 3, when divided by n, it can be neglected, so W=28 (1-6) 


Similarly as above we have 
B(YZ)=E(¥)E(Z)+(4-1684+3487+248°-26B"+8B°-248°)/n , 


and consequently 


ih 1 
me eGR eer 328 +146 '4+248°). (4.3) 


4B7(1-8)* n 





E(B) =B+ 


Now compare (4.1) and (4.3); both of them are divided by 
4 Pe (1-8) - - When 8 approaches O or 1, the values of bias — will 
be very large, though both of the sums of the coefficients of the 8's 
in the parentheses are zero when § approaches O or 1. Figures 1 to 4 
give the shape of the curves of £(8) and bias for different values of 
B . From the figures it is obvious that A is better than B . 

The variances of those estimators are very messy for hand computa- 


tion, and have not been worked out for this thesis. 
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Ve 2OME BASIC ASPECTS OF THE EMA2 MODEL 


The simplest aspect of the EMA2 model is the exponential marginal 
Weoacpmibution of the intervals ix. J ; in point process terminology 
feet. 2} this is the synchronous distribution of intervals and refers 
to the distribution of the interval from an arbitrarily chosen event to 
the next two events. For the Laplace transform of its probability den- 


Peeve function (p.d.f.) Boe (x) , we write 


i (5 )=E "| 
= fe Fee] 6, +2 { - fs CP] (1- @) 6, 


50 Ey 8 Es 45 -8E5 40 ' (5.1) 


tHe 


t-0,)(1=-0,) 


using (1.1). Since the 1.1.d. random variable c. have exponential 
distributions with parameters i , their Laplace transform is A/(A+ts) 


mms (5.1) becomes 








+ A A 
fy (s)= See 1-(,) A (1. -€,) (1-0, ) 
x,‘" xepe Pe" Jefe est 204+ es x x N45 1 
2 <a © 
AB 
This demonstrates that the x. have identical exponential distributions 
as asserted. The parameter A is the number of events per unit time 


@eethe rate of the point process. 


The correlation between Xx. and Xa can be obtained on considering 


the product of x. Prom (lai) with 
CEs 4 WePe Fos 
Kyat?) Es art Es 40 WePe (1-F,)G ye (05 BQ, S14 160,41 52250005) 


CEs les iotisag MePe (1-G)(1-G)- 
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Thus, again using straightforward conditioning arguments, 


BC 9X4 44 )2(BE 5 Es 4 183 *2(BE gE 41 PPE Es 40) Foy (1-6) 

+E(O5E Esa Po OiSaed) P20, (1-8) 

*BOSE Es a4t OP EsEg aot PoE aE nas) Po(t-@, (1-8, 

+B BEE 4y* CCE gaa Cols ata ag) Oo (t-F) (1-8) 

HECBSE Ege * FP E sa aor Po PrE ans Els tba ec) (tb) By 

HEGEL 541 PEE sc Pol sExagt PEs iEsa1E ste 
Ne a ey (1-8,)°6, (1-8,) 

#00G50,Es 41" Pye slate” (i Poses Pier aEsaot ola asEane 
+P ES ae) (AR) Gy (4-8) 

‘(BoE E sas HP EsEs ant PAE, zee OPE r ast Bie, ave a42 


, ae Vane. .5 x2 
*BE ness ag thes asErsa* ibs rc Es voEnag) (I) AR)” 


aieesamplification of this result leads to 


cov(X, ,X : 
Similarly we have 
Urcorr(X, 6X, ,.)=@ (1-8) (1-8) (543) 


Pe | 





By the construction of EMA2, ge fee Loe result for P, (a2) 


equals zero when ome and B,=1 Mae) Bo=1 avi: B,=0 : and will 
approach its maximum value at B, (1-B,)=1/2 - This occurs when 
B,=1/[2 (1-8) ] eee wee Ite fol oe (Ss equals to zero when B=0 ; 
or Bj =1 POY B,=1 ; and will approach its maximum value at B=1/2 
and B,=0 . Therefore, the serial correlations of EMA2 are all non- 
negative and bounded above by 1/4. 


Now the stationarity of the EMA2 process will be discussed. 
Define B(X, )=m(i) (5.4) 


B[X, -m( i) Lx, yom i+k)| =Cov(X, .X. JRO i : eee) . (5.5) 


A stochastic process with a discrete time parameter is said to be 
"stationary" (or stationary in the strict sense) if the distribution of 


cer; 


A mene. eke IS sthe same as the distribution of X, 


e541 ity co aoe 


X ete 6-7 MborievemmGiinite cet of antegers 11,2.....j} 


i+1+k A itjtk 


and for every integer k . This definition is equivalent to requiring 





that the probability measure for the sequence {x, } be the same as that 


of OS itor every integer k . If the first-order moments exist, 


k 


Stationarity implies that E(X;)=E(X; 11) for~all a,kK—-O, +1, +2,...(5.6) 


Since (X57 Xs ) has the same distribution as (xX. ) , existence 


il Ree oak 


of the second-order moments and stationarity imply 
Ogee =O (atkpitatk) . (5...) 
Setting k= -i-l gives 0(1,itj)=O[i-(itj)J]=0(j) . (S28) 


» 


In the normal case properties (5.6) and (5.7) determine that the 


Stochastic process is stationary. 
4 
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A stochastic process is said to be stationary in the wide sense or 
weakly stationary or stationary of second order if the mean function 
premeemescCOVarlance function defined as in (5.4) and (5.5) exist and 
Satisty the relations (5.6) and (5.7); i.e. the mean is a constant, in- 
dependent of time, and the covariance of any two variables depends only 
On their distance apart in time. Obviously, any process which is 
stationary in the strict sense and has finite variance is also stationary 
in the wide sense. In the normal case discussed above stationary in 
the strict sense and in the wide sense are equivalent. 

We have proved that the X, have identical exponential distributions, 


which implies that E(X,) exists and E(X,)=E(X, EOm adel Or amk 2 


k) 
--.- Also we have cov(X x Jet /r~ 
. 4450 ; 


cov(X, 5X54)" ~<a{6, (1-@)-18, (-%) oo 
cov(Ks Fran)” S5[PoltGa) (1-0) 


cov(X, Xs 45)" OF for k=3,4 5scceee 
All these expectations and covariances are independent of time i, so we 
conclude that the EMA2 process is stationary in the wide senSe. 
The independent exponential sequences and EMAl models are the special 
cases of the EMA2 model; these aspects of the EMA2 model are described in 


the following table: 


28, 





Values of G, & G, 
When we set xX, sequence reduces to 
in ENA2 model 


G=63 = 0= 0 EMAL nodel 





Go=bot Oy 0 Po=0s Ey 407Ca44 EMA1 model 


New the adjacent 
elements are indep. 
if keep & 442 nO change 








rts G8, X,"Cs WePe 1 Poisson process (1eied) 
Bo" Of Bet KyPCr44 WePe 1 Poissen process (4-14) 
C_™ OrG= 0 Ky Ss 0 WePe 1 Poissen process (ieied) 


This gives checks on most of the results, for the serial correlations. 


In the 3rd case x. and Xe are independent, so p,=90; Jao X. and 


Jl 


Xo are dependent, so p=B. (1-B.) , which is the same expression of 
Py in EMAl. The serial correlations in the last three cases are all 
memo, Since all of them have i.i.d. elements. 

Also, even the backward model (1.2) could be equally treated to 
produce similar but different results. However, there is no time-reversi- 


bility in the process, in the sense that {XX peoee eX} does not have 


Z 
the same joint probability distribution as {X_)sX_greee Xp} Grea 


Eiaate kK, where k->2. 
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VI. DISTRIBUTION OF SUMS IN ie SEQUENCE OF THE EMA2 MODEL 


In the point process theory of the model, the distribution of the 
sums pes Xt. - +X are very useful; if these can be obtained then 
the distribution of counts, both in the synchronous and asynchronous 
mode, can then be derived. It would, therefore, be a particularly 
attractive feature of the EMA2 model if the distribution of the T could 
be obtained. Unfortunately, it 1S not possible to get a simple expres- 
paeomeor. the Laplace transform of the p.d.f. of T as in EMA]. 

However, a general derivation will now be given. Define wW(s) as 
the Laplace transform of the p.d.f. of the ec. distribution; except 
where otherwise remarked this distribution is exponential of parameter 


X and so wW(s)=A/(At+tS) . Define the triple Laplace transform of the 


me... OL Te ey and es as 


-s,T -s “8 ..€ 
P(s, So 185) =E{o : 3 are! ? rz} for r=4 sc pecccee (661) 


For r=l , we have 


Dy (5, sepsty)=ife 22 22°93 6 ; g ' ¢ 
= 3 ,£°8., 3 + “SC ,°840 
enie 54CoE4 nto "F318 ne ote 2 12 eve. 3 3f(1-8.8, 


safe 1 2s" EF 2b! Bae) (ces (1-6) 


=¥/(084)[0,Ws2)W(s5)#(1- BP, W(E,s, 482) Y(s5) 
*(1-8, (1-8 Y(F,s, 482) Y(6 4455) 


Gal 





Now we relate 9 (S$) 785/83) and 9 (S$) 18578.) using the expression 


r-l 
Tmt they ees O, 
ara ready Mat Ent Es eee (1) P4 , 


Tet BEA Eten Weds (1-8) (1-6). 
Then we have | 


-s,(T ~s +0 jes eof 
rly s8298,)"Bfe 1*\ "red 2 op oor+i seta} e 
tafe tt Fe Aaa Ents 
*Bfo 
= Pinas (84 #Pp54 055) Ws) 6, 


* Prony (4 054 P51 %2)] Ps (1-6) YW(e)4(1-B) (1-8) Yo, 434) ; 


Continuing, we can write 


Dipat (4 54 Sa)™Dyo (54 6PoS4 , P84) Y(sg) oA 


é 
(1-68, 


m8 (Ts * BEG En Emi) “Solty Sate! 


+P, (6, £,84+0,5,78s,)/ 6, (2-8, )is,)+(1-6,) (1-8, YWis,48,)J. 
Dat (249.845 (,84485)= NC fh5. #84) YWR,84 759) C, 


i ran (84 0Fo54 05,465, ) : 


and solve it recursively; the procedure is difficult by hand but could 
possibly be manipulated on a computer. Setting s,=0 and S,=0 , we 


have the Laplace transform of the p.d.f. of TY. : 
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The first few sums have transforms as follows: 


A 
$, (s,0,0)= “ree which implies that Das exponential (A) as is 


expected. For T., we have: 


A. N|O+6,5) At8,8426,5 +6845 2(2-8,)] 
p{e.0.0)= 55 “Caps REOTES eS) ((Gu2) 


os AM26G,8) 
feewe let Bj=1 milion? i reauces to peas A+ 6 mae OF 6 which is 
2 2 


Eeaemmaplace transform of the p.d.f. of T, ieee eee | Ret. J p.8) 


If we let B,=0 , (6.2) reduces to (=)? which implies that T 


is the sum of two independent exponential random variables. 


For T. we have 


2 3 
,(s050)= rs X [(A)(B)4(c) (0)] 


tstG,E+6,5)A+¢,s* 68) T2aG,8 )(A%St6,3 (AtE, 


where A=[s”(265436,6,+670,)tha (50,420, 4° ](1-8,) 





fs )2? 


(6.3) 
Beo( P+ 05+ 05401 Bots” ( Bat pSt3B, Poti Ryt20, eX a (1426, 420, 4)? 


Coss ( G+ Ps* P,05*26 1Co*0re2)#As (b542 f,+2, B)+BaX 
Das? (2p3+# P:05+20,05 #8" (585*6B,0ot ES ACHE EN, 


Tf we let B18, equal to zero or one, we have some interesting results. 


When B,=0 ; (A) (B)=A~ (At) (A+28, 5) * aude Di—o.., so that (6.3) 


reduces to d | A (At2,8) i 


Mo A+B, SS AT Py S 


which is the Laplace transform of the p.d.f. of qT, in EMA] for B=B, 


33 





When 


6,=1, A=(i= 2) (\+ B,s+2s) (A+2 P8 ) : Be(\+B,5 ) (At+2s) At3+ Bs) : 


C= bs (At2s) (At2s+P,s ) and D=(\4+2 ps) (Ats+P,5 e , 


and this will give the same result as above for B=B, : 


When 


Bei, (A)(B)= 0, (c)(D)s+2s+B,s) Ats+f,s)°At2s), 


5 3 : 
(6.3) reduces to i /(Ats)” , indicating that the tx. } sequence are 
isd. exponentials. 


When 


P, = O, A=(i- >) A+0,s ) (A+2 Po) P Bu\(+f,s) (Aist 6,8) 4 


c=B (ats) Ats+ fis) and po(AtG,s)" +268) 


So that (6.3) reduces to 


*2 0-8 
Ce 
Gas) Cap pyOrstts) 


which means x, and X., form an EMA1] model, X., is exponential (A) 


and independent of x) and X, “ 
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Vit tie vOLNT DISERIBUTION OF x. AND X, IN EMA2 


it] 


We now discuss the joint distribution of Xx. and Xa which will 
be a bivariate exponential distribution. Several authors have discussed 
bivariate exponential distributions, including Downton (1970), who makes 
some comparisons with those of Gumbel, Moran and Marshall-Olkin. The 
distribution to be discussed here does not appear to be one of the 
earlier ones, although it is fair to say that in common with earlier ones, 
it is not the 'perfect' bivariate exponential. 

The double Laplace transform of the joint p.d.f. of xX. and Xo 1s 


I 


easily calculated using (1.1); the required expectation is 


ia 
eel 


E fois Saha] of enn 
“E fe ap a ro 2Pa+) Pan (o 88 465~Sa( Pog ag PE 172! 8.0, (1-2) 
tefe Roe 1Ea“Salfotan “ie vsa%a9) B(1-F) (4-8) 

+E fo (Coes PEs as)“ PeSaf a 6,0,(1-8,) 

+5 fo’ (Paes Preset) Pa' Pot set CiEsaa’ a2 2 

+5 foo rai" Prenat Pa CEs the usaitass) B, (1-B,)°(1-8,) 
A aaa AC NICE CD 

‘a Robs tess *Este) wealbafaas* iFssa) 6, (1-0)°(1-B,) 


+E 


45 fo (CoE ,* PE 541? a0) “sol PoEs ea Pye 442°C443 ‘ (1- B,)? (is 6, © 
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which can be written 


.- Xs 44 Ps ra WR 21) Bo WC 85) 


Xye 


+0,0p(1-B Yb 59)W(Gs9)*W(R,s4+6,8,)] 

-+B,(1-8) (1-8) [WGs,)W(6,s,)We)* WG,8,+6,) Ws,)] 
+6r(1~6,)°[W( B24, (B,s,)] 

+B (1-6,)°(1-6 LW 0,448,855) V( O25) V(55)+Y/(8,5,+8,55) W(s,+0,5,)] 
+(1-6,)?(1-8)?7[W( Pe W Css" Bate) Vler*Fie2) Ve.) } 


_X{O+8,94#8.0,) O89, 46,99) O49,)#8,(1-8,)9100[BN14G)-G,A-B5,-0, (1-8.)59]) 
+0484) \A*S, (A*8 5) Tie, F9,8,) ATS Fen 
(7.1) 


We note that (7.1) iS not symmetrical in Sy and S, poand this’ 1S Lo 
be expected since the process is not time reversible; this is one 
feature which distinguishes it from earlier bivariate exponentials. The 
backward moving average model (1.2) corresponding to (1.1) has the joint 
iiterval distribution which is specified by (7.1) with Sy and S. in- 
terchanged. 

eeee wo licie form of the J7Oint distribution (7.1) can be obtained 


directly, rather than by inversion of the transform which is less infor- 


mative. By the structure of the model the joint distribution of 
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(XX. 14) Ha WieteMeeron Ene JOInce distributions of 


(BE, PoEver)» (ob PoEres* Este)» (Eso Poss Eso tby ag) » 
(BES*CEs ere Esa)» COE TEs ere Err Er sets ag)» 
(BEs*+PEx ay Pols sat Ener)» (RExt AE argos eset Ege)» 


(GEC Es a tensor Cobses)® (CEs tls as Eves Coban AEs 12*Fy 29) 


With corresponding probabilities 
oe B, 0, (1-05) « 0p (1-0) (i-G,)s (8, (1-8)» @ (1 B,)° (1-8), 
B16)", A (t-0,)°(1- 8) + BUI BMt-G)s (1B) 8) 


These joint p.d.f.'s can be listed in an obvious notation as follows: 


Eases vPotaag Mba)” exp(-x/B2) exr(-ry/pp) (xpy>0) 
= 2 = 
AA eed MCa6G~bo)) 
-oxp(~rx/¢..)[exp(-1x/p,)-exn(-rv/p..)] (x,y >0) 
| 2 » . pec 8 
fetta» Colaay ee © ba {exp[- oe Ay/b.)]exp(-Ay/n)} (6x>6,y0) 


= \- [(4. fy] 74 
* be Cotas thEgaotsag A (1 )(B, 6) 


vexp(-Ax/Pp)} oxp(—y/6,)~exp(—\y/6, )~exp(~Ay)texp[-Ay (1+ 3, “a (xey>C) 


27] 





Meeheritemis are mone difficult. For example, take 
(CEs Egat CoEs rit GEs 40) 
Let x (Fy*0Enege TOEs at A esses Gaye tes 
E440" “/B4s Exa2ly-2)/ Cos € a= (8yy-G2)/ Bol/Bos 
and the Jacobian equals to | 


2 fe 2 pai 
1/856,» and BQx=G5E,+0PiEs24°0oC, th, y~0,2 which inplics that 


Bx> Byy- hz => Bz > CY~EX =p 2? ¥~x/64¢ 


Thus when Box>By , we integrate with respect to z from zero to y , 


but when 8 x<B LY , we integrate with respect to 2z from Y-B.x/B, ee. 


2 


y . Hence we have: 


when Box> Gy 70 


2792 2y=1 
f 3 = Ay reA =) ee ® 
tatatOrEsnssCrPaat* esac” FOX (Bo- Gi PatBi)rexp\2/69) 


{exp [-Ay(6o~0, )p5 Foxe (-Av/0, ) 
when ysGx >0 
= CEs ate PoEsn1* ee DX Ooh 00) exp(~r¥/B,)° 
f exp [-Ax(1-F/6, \/6,pexp(-Ax/6,)} 
a (BoEstC Esa? Cobar trEssotes+3? 
Let eB E +B Ess AEs atte gscteaage BM RExsar MExagt HH 
Eqn gts Cy40" 2/Oy1 Eas ret)/ fas Ex=fe-( Oy B2- Py) BI/ Pos 


2 2 - 
and Jacobian =1/85 044 also Gx PoE, Py PoEs 4s Pabst yr~ Pie AM which 


he. 





implies that (x > Gy~G2-0,¥ >Pu > By-P, 2-0, Sw >y-2-Bx/B, 3 for 
y<2thx/ G,» integrate w from zero to y, for y> at Bx/ 6,» integrate w 
fron y=z~ Bx/ 6, to ys; and in the 2nd step, since y< at Ax/e, implies 
that z>y-B.x/6, ae hiy< PoX. integrate % from zero to y, if Cy > Box, 
integrate 2 fron y~Cx/B4 to y3 also y>a*hx/6,>2 <n FX/ By» thus 4f 


Cx<b,ys integrate z from zero to y~Px/p, » if Cx <x f(z)=03 hence 
(x9y)s 


for the expression of 


£ 
PoEs* Esai? Coes es iE 42"Fs13 


(xay)=(BA)* [(B2+67-P io) Po*Cy-Po)] exe (*/8) 
fer [ay 1/p_° by /p pene oly) ~<xo(Av/6,) +02 [-ar(t41/8,+ 0/0340) 
wen s>veBir/pys yy >0,x>0 | 
atiav (BA)? ((6358%-6:8,)(03"G,-Ba)] Hf exolarz/pp)-oxn orxCt-6/6,)/0,]} 
{exm[-ay(141/0 *f,/p5-1/82)]~exp(-Av/6,)} 
when u<y-Ohx/pit Byy¥< 0x 


£(xsy)= 0 
when zcy-Bx/ays By > P70 


e(xsy) DX [(8240-G.)(1~0,))™ 
fexp [-ax(1-0,)/B, Jexp (ay) exp [(-64) etx e/6,)} 
~(AB,)” [Bo#85-P1B2) (B2*Pr-Po “erLAcee 2) 
~exp [-ax(1-8,/8,)/B,) ep [Ay (141/08 + 6 /5-1/8,)] 
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The rest can be derived in a Similar way. We thus see that the joint 


pea.c. of XX) will be continuous in both variables but will have 


different analytical expressions over the regions Bx°Biy and 
B.x<Biy ; there appears to be no compact analytical form for 


f 6 [a eee ats lS UnrOortunate because 1t makes it difficult to 
1 itl 
derive maximum likelihood estimates of the parameters i and 8 in 


the model. 

Different bivariate exponentials also can be compared through their 
conditional properties and so we will derive these for the present con- 
tribution. Conditional p.d.f.'s are not succinct enough, and so we 
concentrate on conditional moments. These may be obtained from (7.1). 


Remeeinstance, to obtain Ce) prcitecrenerate with respect to 


Si set Salas mitliety oy ==), anvert with respect to S. and then 


divide by the marginal (exponential) density of Xa . Thus 


B(X4/X, 445 t)= A {146 B+) -6, (iB, )exn[-rt (10, )/6,)/(6,-B,) 
(6-2, +F,B.dexo[-Av(1~A,)/0,1/(6,-,)] 


Examining this regression function more closely we see that E(X, |X, =t) 


i+] 
-l 


‘is equal to 2 fens B.=0 and By equals either 0 or 1; otherwise it 
-1 
increases exponentially from (B,B.+B,)A to the constant value 
-1 : 
(1+8,8.+B,)A as t increases. But when Bo=1 and B,=0 ; 


B(X, |X, )=t)=3/2 Vireciersetic Mdximum value for large Ct. 


The conditional moment E(X, |X, =t) can be obtained similarly by 


i-l 


interchanging Sy and S,5 - 
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ViTI. SOME BASIC ASPECTS OF THE EMAk MODEL 


By the constructions of EMAl and EMA2 model, we can write the 


general form of EMAk as: 


xyes MaPs Py 
“Hest Fen 5 44 : rape (1-6 B24 (8, 1) 
“EtG. e141 Oenk 142 ; Wop, (1m G) (i- = )Pr2 
“B8,+0 Esai so tes Die ree One WePe (1-6) ( See eee (1-0, ) 


Methods of mathematical induction will he used to prove some basic 
properties of the EMAk model. 


meeethne general closed form of EMAk (k=1,2,3,...) is 


is J 
et] an 
x= = 8, .E4,, TT ieee (8.2) 
agi YY et 
(k+1) ' 
where Bo and I. are defined to be identically 1 for all i; 
ae 1S an 1.1.d. sequence of Bernoulli random variables 
; (m)_ : 
with I. 1h ior (1-6), O otherwise for all nm; 


Petomene Serial mumoer of the Pee lenen t Gi the Serics: 


k is the order of the process; j and n are indices. 
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Proof: When k=l 


x RE? A orth ; 
alas wep. (, 
ee wep. 1-8, 
When k=2 


Bye) COyEe)-() 
1-R2,117 +0 s+it aes ‘r * (Es aatg Tg Ly 
=(3€, , WePe &% 
“E+E ass HePe (1H (3)0, 
“BE tesa tegaoe WePe (1-()(1-8) 
Assume the result is also true when k=m then, when k=m+l 
. s- gfeFlon) (¥en). 
r alti _ \ercen 
“= a eS ats = +0 tle S445 Tay 


m1 


7 2 & Paring’ ” Tr T x(etieten), 


This completes the proof. 
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2. The distribution of the intervals {x, } are also exponential. 


2% “BA, 
Proof: fy (s)#E(e Ne, lets 
di 


% “she 0, 2.75, 0 
Barat af: 196 +2(0 eva h : Aug AB) toee. 


+0( ve ea ore Na Fi.) (AB ove (1-8, 


ry, XA+B8-4 X77 (1-8, owe (t-B,.) os 

“MS * OBS Oe as) * DRTOF TOT _ 
When k-1, 

Pca X(1-6,) 


fy, (6)* 555.5 6,8 * OsOd 68) 7 nee Soyo ngage am 


When k=2 
by ro, Xm “0d +60) X Naf.) ) 
a pal bot TORS)” waa" ea J xh 


When k=m-1, the last term of (8.3) is 


0 (4+6.) (1-6..4) ece (1-B,) 
Crs ORs) OFF 78) 00 O88 





Assume the result is also true when k=m-1, then, when k=m the last two 


terms become 


_Xa-BG- Graadoee Bo) Nei (1-B.)(1-Bog) vee (1-8,) (8,01) 
OF P,E) OB Sone OBS) ARS) OMB COB, go) 000 ON BLS) 
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but all the terms before these two are still the same as k=m-], thus 


Simplifying (8.4) gives 


(1-8) (1-4) (1-B)[ 8, Ovts)#C1~8,)] 
~ (ts) (\*6.,8 OTB 48) oe e CN 658) (K*Gy8) 


Nf, JB) eee (1-f,,) 





25 
which 1s exactly the last term of f. (s) when k=m-l. Hence, we 
A 
proved that if the result is true when k=m-1, then the result is also 


true when k=m. This completes the proof. 


3. The jth order serial correlation of EMAk is 
1 2 2 5Glc 
ot TL ir jo 0 for 1i%j3&k 
J ie = Mitt mi (1 ge ee “( ktien ), 
ic 
= .0 for 3 


where Gy=i rand P41" 06 


Proof: By definition 


cov[x¢) .x(h) 
Pea aT 
re 


_ lef rf6 fo) fe) 
fran vex [X pj van [x Pali : 





(k) 


where xX. '"'s are intervals of EMAk process, and have been proved to be 


marginally exponentially distributed with parameter A . Thus 


arava) Pals = 


o(_f (ky (kK) 1 (x*)] p Cs) 
oH) x {3 xt ite Ses ve : 
Since a and aie are probabilistic linear combinations of i.i.d. 


exponential (A) random variables ¢€. and €... , and 
it it} 
B(E Ey 4,)°B(E, ECE, ga )et/ 0s 
icat5 " itj 


the only non-zero term of ps will be the sum of 


BY [BCE 34 23)7B(E, £5, +3 ) =BY(2-1)\“=B, 


where B is a combination of B. and (1-8. ) Pe Oi 7 7S eve 


Hence when k=l, j=l, 


, 4 4 <4 \e are 7 
ee 4 es 4 | Er a = 
= Rats i (4 Catan itl<-j=<i up Patan 


=6, (1-8 hy (1-8, =O, (1-8) 


When k=2, j=l, 


7). 36.7 TT Se Oa Tr (1 en) ss 


=((1-05) (4-0, )#0, (1-03) (1-8) Py (1-,) (1-8,) 
~ BQ (8-Q3)+By (4-8) C1-G,)9 (1-0)=[ (2-09) 
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When k=2, j=2, 
- 1 
C5 PAS TL 1-F-nd@inn TG an) 


=B,(1-G,)Qp (1-6) (1-8 = (4-8) (AG) 


When k=h, assume the result is also true, then, when k=h+tl, j#h+l, 


(Ik). hejtt itje 
— = Cnt “4 0 5 ee Oe TT e “Cheten) 


fof itje4 
+P sete(ntteset) TT t-Gritetam Cnet tinged qT (1-6 steten? 


Ietlej 4-4. itjel 


a oGatut 1) (i- Catton) begus 4 TT (1G igen) 


This completes the proof. 


All the correlations are non-negative and bounded above by 1/4. 
Proof: From above 


itj-1 | 
0. SPs T (1-Gortom) Pitt ajo ue (1-Gpten) aoe 


a (1-4) e068 (1-8 423) 
4 ki (1-8 Ga gag (toR) (1-8, 4) ees (1-85) 


+B o(1-B 1B Gen gant 1-Beag dows (1-G, 54) 


#0,44(1-0,.)« 00(1-B. 40) G 1G) (1-004) eve (1982) 
+0,(1-G.)(1-G dove (1-Perg)(t-B DG ag dee © (1-8, )» 
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1 
When k=l, j=l, 0 =B, (1-B,) r 


Min. value=0 at B,=0 or 1, 


max. value=1/4 at Bi=1/2. 


(2) 


When k=2, j=l, py ?=8, (1-8,)-(8, (1-8,)17, 


Il 


Mon. Value—O0 at B,=0 or Bo=1, 
max. value=1/4 at B, (1-8) =1/2. 
Rigen) k=2, j=2 p\?) = (eas Cass) 
‘A : 2 Z it Dre 
Mune valuec—O at B.=0 or Bi=1 Cr Bo=1, 
max. value=1/4 at Bo=1/2 and B,=0. 
Wren 1<j<k, Min. value=0 at BOs Ms pte ly eee pk 
or Be Mak,Or kK-],...0r k=j+1, 
Max. value=1/4 at Ble and B =O iene irae 


Wihewe weM=k phi, Ka 27 .<.0 7 ] s 


peeeDetine the Heke f element of EMAk to be 


CEsate Were Qi, 
ove i~ 
Cesar tates se! Wore (I~ Fh) Fheat 
(k) 
Mit C Ear tente stat C2543 weve (1B) Go) Fiat 
CEs sitG sEgaatenet hls ancet PEs Earnest WePeAnG) AG ag? 
oo0(1-G,) (1-8; ) 
Geemcerine the 21th tlement of the k+lst order process to be elie ; 


then we can write 


(iH), 


: (0B 414 10,41 st2ye0e) (845) 


(k) a 
Geries* Xg44° os a 


—_ WePe Fras 
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Proof: 


when ket, xU2eBE Wore @, 
CEs tesae ° WePe Inf 

Then x86, , WePe F 
“HEC Ess» WePe (1-0 )B, 
“BE AOE 1tEs ios wpe (1-8,)(1-8,) 


(2) 
When ke=2 $ Kea anaes 9 


- DE, a tREs sos 


WePe @, 


Webe (1- 62, 


OE ut AEs otE yas WaPo (1-0) (1-2,) 


Then 


x nBE, WePe e 


“REE » 
“QE + BE ait ereee 


3 BE5* CEs 4" GE 112*e443° 


WePe (1- b> 
Wee (1-@) 1-8) P, 
WePe (1-8) (1-@,) 1-04) 


When k=m-1, assume it is also true, then when k=m, do the same job, 


Pemeeget exact the correct result, this completes the proof. 


Note that this expression is not convenient for the purpose of 


examining the properties of the EMAk proce 


dependent; however, it may be used to gene 
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IX. CONCLUSIONS 


Both estimators of 8B are not very good, since 
very large when 8 approaches to zero or one. 
nice when pelomenecne interval (Ol, 0.8)x 
Estimation of 8B in the EMA2 process is rather 
it is impossible to get the unique value of the 


In successive stages of queueing lines, all the 


the bias terms are 


But 8 looks pretty 


dtiiacule,. because 
estimators. 


waiting time “(waiting 


time in the queue plus the service time) in each stage will not be 


independent; this is the basic purpose of constructing this model, 


the size of the order k depends on the number of stages. 


The general expression of the EMAk model is not 


convenient Lor the 


purpose of examining the properties of the EMAk process, since 


«09 
alk 


the ix, } sequences. 
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and c. are dependent; however, it may be used to generate 





APPENDIX A 


METHODS OF GETTING JOINT EXPECTATIONS 


The standard way to calculate the expectation of two or more jointly 
distributed random variables is to integrate the function with res- 


Meeteto the joint p.d.f£. of the random variables; e.g. 
er)= Jy Sy oy fy y(xuy) axdy. 


This is not convenient for the expectations we require. 
In the EMAl model, a better method of getting joint expectations is 
to write out the expressions from the basic construction and compute 


; . 9 
them directly. For example: X,=00, Wepe G 


GE tls 44 Wepe (1-8) 
2 pert 
Kear Ert1 wave Pi 


eee 2 . o (1-(S)> 
“OEn as eres uesto tase NeP (1-8) 
nus 


2 2 ere of 
Kyxr sO E seat eee 
3 2 ot e 3 WeDe B(1-0). 
H&E ses Ci +t 
2 " 
“BE Ear tee gateaag PE seize Pe oe 


a 2 2 Zug 2 Ae te ea 
OE Enns t2O Ese aes ao PE pas Ey Pyro ese ii2 C5 410442 
Wepe (1-6) 
By direct computation, we have: 


BK EX ag) H/, ss MePe Bs 
=(120°+200 '4288°)/p? Webs 6° (1-6), 
m(120°4568'#760745207H48)/>> were (1-8), 


=(4p 4288448 B42 04160" -H4) / 2 WeDo (1p)? R 
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Combine the above gives 


BO KS a4 Ef ag) Catpr20e? -20p749¢°=126")/d. 


All the expressions listed in APPENDIX B were computed in this way. 
Mawe the derivative of the Laplace transform of the joint p.d.f. with 
respect to Sir and then setting S. equal to zero will also give 
the joint expectations, e.g. Lawrance and Lewis gave the general ex- 
Meessi0n Of Laplace transform of the joint p.d.f. of r adjacent 
matervals. [Ref. 1, p.17] 

Converting it gives: 


EK 


fy x 


(S48, 9S458),)™ 
4 settasokgag Pe" 304 


be (\4+6s8 1788. )(At6s,, p88.) At88.*88),) 
(ips 1) Os 2) O07 is JOP B52 Oe," 683) (NS 4765) (X48), 






Take the derivative of this with respect to S) twice, S. once, 


S, twice, and Sy once. Then set S,=0 (iat 2,73,4) to get 
2 e s ° 
B(x? > ee .« Dee ) . Note that when the order of the derivative is 
1 it] it2 i+3 


odd, one should change the sign of the expression. This is a messy 
job by hand but one done easily by computer. 

An alternative way is to use "cumulants" or “semi-invariants" 
(bets. co, D.255 and ll, p.55-93]. Let L be the Laplace transform 


Piece FOlnt p.d.f. and L*=log L. Let L* 412 denote the derivative 


twice, s once, Ss once, and s&s Ewiece. 


— 
of L* with respect to s > 3 4 


iL 
Z Z 
Since L*=L,/L and L* ={L,L4L,) }) /L and L(0O)=l1 , these imply that 


 L*, (O)=L, (0)= -E(X) and L*(0)=L, (0)-{L, (0) 1 “= var (X). 





we Gece from this 


If we denote ee mao) = and FE J ym = 
Km ei) Es 


relationship the following: 


+ Ko ; 


E 43 


= kK 


Lae ae tt 


hee 
Bog? Kok y (2, 41K) Ky 


E. 4" “Ky onK, (2K, 47 “Ks 


i Koo) 


2 Fa ie x y by 
Boot Kn ot2Ky (Kyo ato AREY Ko Kg tA (4 Koy 4 


: . tb 
By 3° Ky oy (Kyat Ky 943K, (Ky ot Kr (Ky 1X go) #4 5 


Bai" Kg thy (Kggt Koy 43K y Kot By (Ky HK bog) Ky coe 


where = oe cS os SS « G3 en}! 
Ry Kym Kogt ~t/) = -Ey™ “Eig, 


Also K, 4=(60°)/>%, 
Rag Kooi ln 
Ky o%(26°~28)/»?. 
Ko (2-28) /?, 
Kyot(68--68")/\, 


Ko3™% 30" ~2/)), 


3) Pe 





APPENDIX B 


Polson USeEUL JOM EXPECTATIONS 


- + i Zz 2. 
E(X.X..4) (1+8-87)/ 
2 = = fa? 3 ee 
E(X2Xs 44) (2+48-28 Bea, 
2 )=(24+2B-2B°)/d°- 
E(X.X. 1) (2+28 8 Me 


E(X2X2,.)=(44884887-148°-28")/ 2". 
3 = yoo (<a & hy 

E(X;X. 4) (6+24B-188*-68')/d°. 

3 s Eris Fu i 

E(X Xe) (64128 6B )/d 


ECKEXS ) = (12436846087-488°-368"-128°)/2°. 
E(x: ye 2 ye(izsz4ee24e?# 126% -403"-128°)/2", 
E(X2K2,, )=(364108841808"421689-3248"-144B°-968°)/ 26, 


E(x: x. 5 44968-2482-2482-48R"4+248°-24B° )/d”. 
244+248-2482+48B" -48B°)/A°. 


1 0046006-12082-12082-1208"-1208°-1208")/2°. 


onan 
E(x. x. aes ( 
E(x, 5. ( 
vine. 


E(X. Xe a 1 Dea Skee | Oe yeas 


E(XEX: \=(484+9684+9682 +488 °4488"-2408°-48B° )/A°. 


a 
e(xtx2,, )=(4sv19094320%-24092-19284-14482-4888)/26 
Flysy4% ye . 
(XiX ey) (5764+2304845184824+864082+120968"-167048°-80648° 
SU? POS VON 





= 23 3 
E(X.X. 15 %54,) (1+B-B°)/A°. 


E(X2X i 41 Xe po) = (2468-48°-28°)/2?. 
E(X. X24.) Xi 49) = (2468-887 428") /2". 


Z = - Pg 3 3 
E(X Xs 44X55) H(244B-28%-28°)/d°. 


me 2 S34 QR QH_ 5 6 5 
E(x? aXe AX +9) (441684368 968749687 -88B°+368°)/dA°. 


E(x 2X 4X3 42) (4412B8-487-148°4+128"-6B8°)/2° 
E(x. Xe ae X25) = (444842087 -~208°+88°-128%)/d°. 


E(XEXS X2 .5)=(8+328+4887-568°-408" 488° +88") /d”. 


a | 


E(X.X5 4X. 45) = (64684428%-548°4+168°-108")/2d°. 


E(XeX: x =(124+968-3687+608*-1568%+2048°-1688°)/2§&. 


542! 
XS 5) =(1241688-2768°4348R°-3128"4728°)/2°. 


1 | 


E(x, re 4] 


EO, ; (are Sane an 3)=(44+1684+68B*-1328°+1048%-160B 541248 °+4B7 


SQUIB? yee 
7 
BAS, 


. | | 
Xe go %e gg) = (4420841 287-288 °-248"- 48° 4568-4487 


+128°)/2°. 


| 
AKERS Xa oks )7(44248-44874808°-1088 "46085-2085 +1287 


-4p°)/d°. 


E{X. Fare p4o% sag) = (44168448? 36874528" -648°+368°- 168” 
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2 e een Bho 4 5 9Q6 5 
E(XS XS kok yas) (2+88-48°4+28°-148'+10B8°-2B°)/A°. 


2 = 2 3 boas _ 7 8 S 
541% 442% 443? (2+88+28°-268°+268'-6B8°-12B8'+8B°)/A°. 


2 = 922. 3_9p4_925 5 
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4 = ae a S 6 6 
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